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Spontaneous decay process of an excited atom placed inside or outside (near the surface) a carbon 
nanotube is analyzed. Calculations have been performed for various achiral nanotubes. The effect 
of the nanotube surface has been demonstrated to dramatically increase the atomic spontaneous 
decay rate - by 6 to 7 orders of magnitude compared with that of the same atom in vacuum. Such 
an increase is associated with the nonradiative decay via surface excitations in the nanotube. 

' Theoretical prediction of the Purcell effect in 1946 |l]] has stimulated series of works, both theoretical and experi- 
mental, aimed at the detailed investigation of the phenomenon (see recent papers J^, [| ||, ^, |) and references therein). 
The effect lies in the fact that the spontaneous decay rate of an excited atom essentially depends on whether or not 
the atom is located near optical inhomogeneitics and interfaces of media with differing optical properties. Depending 
on the specific configuration of inhomogeneities (interfaces), the atomic spontaneous decay rate may both increase and 
decrease compared with that of the same atom in free space. The Purcell effect took on special significance recently 
in view of rapid progress in physics of low-dimensional nanostructures. It was shown to be of great importance for 
microcavities Q , optical fibers || , photonic crystals Q , semiconductor quantum dots . 

InRef. §, it was first suggested on the basis of the model of an ideally conducting cylinder that spontaneous decay 

S process of excited atomic states near a carbon nanotube (CN) might possess nontrivial peculiarities. However, the 
ideally conducting cylinder is not quite an adequate model to describe the optic properties of real CNs. It is our 
purpose in the present paper to give consistent consideration to spontaneous decay processes of the excited atom in 
the vicinity of CNs. We calculate atomic spontaneous decay rate variation for infinitely long achiral single-wall CNs 
of different radii and demonstrate that the decay rate may dramatically increase due to the nonradiative decay via 
CN surface excitations. 

Quantum theory of the spontaneous decay of excited atomic states in the vicinity of CN requires the solution of two 
fundamental problems. They are (i) the problem of the macroscopic description of the optical properties of solitary 
O , CN and (ii) the quantization problem of an electromagnetic field in the presence of CN. We solve problem (i) based 
upon the model described in Refs. || |{J. According to that model, CN is changed by the infinitely thin anisotropically 
y—i \ conducting cylinder with effective boundary conditions imposed in such a way that the field at a distance from the 
J> ■ cylinder surface be identical to the actual electromagnetic field excited in the system. In so doing, only the axial 
conductivity of CN is taken into account and the transverse conductivity is neglected for the following reason. The 
axial conductivity forms by intraband and direct interband transitions of 7r-electrons in CN, while the transverse 
one only forms by indirect interband transitions p!l| . The intraband transitions dominate at lower frequencies 
and the interband ones start essentially contributing to the total CN conductivity at higher frequencies. However, 
■ the contribution of the indirect interband transitions, being stron gly suppressed by depolarization fields, is always 
smaller compared with that of the direct interband transitions |Io[ We use the dispersion law for 7r-electrons in 
the tight-binding approximation with allowance made for the azimuthal momentum quantization. Energy dissipation 
is taken into account within the relaxation time approximation. The spatial dispersion of the CN conductivity is 
neglected (see Refs. ||, |j| for its role in CNs). 

The quantization procedure of the electromagnetic field, problem (ii), faces difficulties similar to those in quantum 
optics of 3D Kramers-Kronig media where the canonical quantization scheme commonly used does not work since, 
because of absorption, the respective operator Maxwell equations become non-Hermitian. A standard approach 
£j \ overcoming these difficulties involves a noise current term incorporated into the operator Maxwell equations p"2] ] . We 
L" ' use the analogous approach to quantize the electromagnetic field in the presence of CN. In this case, the noise current 
. £^ \ becomes the surface one and, therefore, may be incorporated into the boundary conditions for the Maxwell equations 
rather than into the Maxwell equations themselves. As this takes place, the effective boundary conditions for electric 
?H ' field and magnetic field operators in frequency-domain space take the form 
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where r is the radial spatial coordinate, R cn the CN radius, n and e z are the unit vectors along the external normal 
to the CN surface and along the CN axis, respectively, a zz (io) is the axial dynamical conductivity of CN (see Eq.(36) 
in Ref. B), is the axial noise current operator. The latter one is expressed in terms of 2D scalar bosonic field 
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operator /(R) (Re CN surface) as = ^ tuvRe a zz (lj) / n /(R) and is responsible for correct commutation relations 
of E and H operators E3] . Homogeneous Maxwell equations along with boundary conditions (|l|) describe quantum 
electrodynamics of CNs. 

Let the excited atom with an electric dipole transition allowed be located in the vicinity of CN and let the atomic 
dipole moment be oriented along the CN axis. Following the quantization scheme above, one can by analogy with 
Ref. [D obtain the Volterra integral equation for atomic decay dynamics from the upper stationary state to the lower 
one 

C u {t) = l+ f K(t-t')C u (t')dt' (2) 
Jo 

with C u being the occupation probability amplitude of the upper state and the kernel given by 

K{T) = t3t^L r du **G„{t a ,t a u,) r_ i(w _ Wa)t 1 

h J 1{LU-LU A ) L J 

where \x z and tu A are the matrix element and the frequency of the atomic dipole transition, respectively, r A is the radius- 
vector of the atomic position, k A — lo a /c, and G zz (r, r A , ui) is the axial component of the classical electromagnetic 
field Green tensor in the presence of CN. 

For further calculations we need Green tensor components G az (a — r, ip, z are the cylindrical coordinates associated 
with CN). We use the representation G az = (k~ 2 d a d z + 5 az ) g with g(r, r A ,ui) being the scalar Green function of the 
electromagnetic field in the presence of CN. The latter one is expanded over cylindrical waves in terms of boundary 
conditions (|IJ), so that at r,r A > R cn one obtains 

_ Rcn_ ^ ipv f v 2 I 2 (vR cn )K p (vr A )K p (vr) %h 
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where go= exp(ik\r — r A \)/4Tr\r — r A \ is the free space Green function, (3(uj) = 4ni a zz (cu)/uj, v= \fh 2 — k 2 , I P (X) and 
K p (X) are the modified cylindrical Bessel functions. Integration contour C goes along the real axis of the complex 
plane and envelopes branch points ±fc from below and from above, respectively. For r,r A <R cn , Eq.(^) is modified 
by the replacement I p <-> K p in the numerator of the integrand. 

Consider decay process in the Markovian approximation. Then, factor [exp(— ixr) — l]/ix in Eq.(^|) is changed 
by tt5(x) + iV(l/x) (V denotes the principal value), and one arrives at the exponential decay model with K(r) — 
— r/2 + i Sou, where T and 8ui are the decay rate and the Lamb shift of the upper atomic level, respectively. For the 
outward atomic position (r A > R cn ), the decay rate is written in view of Eqs.(p|)-(^) as 

— = £( 1 = 1+ V [ PAV 4 A ^(vARcu)K 2 (v A r A ) 
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with (3 A = [3(uj a ) , v A — y/h 2 — k\, Y = 4fc \ \ /j, z \ 2 /3h being the free space decay rate and £(lu a ) representing the 
influence of CN. For the inward position (r A <R cn ), Eq. ((HJ) is modified by the simple replacement r A <-> R cn in the 
numerator of the integrand. Note the divergence of the integral in Eq.(|^) at r A =R cn , i- c. when the atom is located 
directly on the CN surface. This divergence originates from the averaging procedure over physically infinitely small 
volume when describing CN optical properties. Such an averaging does not assume any additional atoms on the CN 
surface, to take them into consideration the procedure must be modified. Thus, the applicability domain of our model 
is restricted by the condition | r A — R cn \ > a, where a = 1.42 A is the interatomic distance in CN. 

The decay of the excited atom interacting with medium may proceed both via real photon emission (radiative 
decay) and via virtual photon emission with subsequent medium quasiparticle excitation (nonradiative decay). Both 
of these decay channels are present in atomic spontaneous decay rate V described by Eq. (^) . The problem of the total T 
partition to radiative and nonradiative contributions is not trivial. Radiative contribution T r was estimated by using a 
Poynting vector approach for the atom near a microsphere in g] and for the atom inside an optic fiber in || . Following 
this approach, we estimate the spontaneous radiation intensity distribution I(r,t) = \ F(r,r A ,uj A ) | 2 exp(— Tt), where 

F a (r,r A ,uj A ) = -4mk A Li z < G az (r,r A ,uj A ) V ■ duo \ , 
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at large distances |r| — > oo. In so doing, the second term does not contribute and the contribution of the first term is 
easily found by the stationary phase method (Til to give in the spherical coordinates (|r|, (j), 9) associated with the 
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The relative contribution of the radiative decay is then given by 
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Figure [l] shows the results of the numerical calculations of factor £(o>a) according to Eq.(||) for metallic and 
semiconducting CNs of (m, 0) type ("zigzag"). The atom is located on the symmetry axis inside CN. The presence 
of CN is seen to drastically accelerate spontaneous decay process of an excited atomic state. Frequency range 
0.305 < huj A /2^o < 0.574 (70 = 2.7 eV is the carbon overlap integral) corresponds to visible light. Lower frequencies 
huj A /2j n < 0.305 correspond to infra-red waves emitted by highly excited Rydberg atomic states. At these frequencies, 
the large difference (3 — 4 orders of magnitude) is seen in the value of £(uja) for metallic (m = 3q, q = 1,2,...) 
and semiconducting (m ^ 3q) CNs. The difference is caused by the Drude-type conductivity (intraband electronic 
transitions) dominating in this region, whose relative contribution to the total CN conductivity is larger in metallic 
CNs than in semiconducting ones |J |To[ O. As the frequency increases, interband electronic transitions start 
manifesting themselves and function £(u) A ) becomes irregular. At high frequencies, there is no significant difference 
between metallic and semiconducting CNs of close radii. Function £(u> A ) has dips when uj a equals the interband 
transition frequencies; in particular, there is a dip at h,oj A = 2^ for all CNs considered. It is essential that £(ui A )^>l 
throughout the entire frequency range considered. This lets us formulate the central result of the present paper: the 
spontaneous decay probability of the atom in the vicinity of CN is larger by a few orders of magnitude than that of 
the same atom in free space. In other words, the Purcell effect is extraordinarily strong in CNs. This is physically 
explained by the photon vacuum renormalization: the density of photonic states (and, as a consequence, the atomic 
decay rate) near CN effectively increases (p e ^(ui) du= £(w) ui 2 dui/irc 3 ) since, along with ordinary free photons, there 
appear the photonic states coupled with CN electronic quasiparticle excitations. 

In Refs. H ||, the possibility was shown of the existence of slow surface electromagnetic waves in CNs. In Rcf. [§, 
such waves were shown to be responsible for the strong Purcell effect for the atom in the spherical microcavity. The 
results of the present paper agree qualitatively with those obtained in Ref. Q. Quantitatively, ^ ma£ c~ 10 6 — 10 7 for 
CNs while £, ma x ^10 4 for the microcavity i. e. the Purcell effect is much stronger in CNs. It is worth noting that 
for the atoms with large enough Tq, there is the risk of going beyond the applicability limits of the two- level model 
and Markovian approximation. Then, in the first case, the mutual overlap is possible of the levels due to their strong 
broadening and, in the second one, the atomic spontaneous decay may become nonexponential so that the problem 
may require the numerical solution of integral equation (||) with kernel (||). 

Figure]^ shows £(ll) A ) for the atom located outside CN at different distances from the CN surface. The qualitative 
behavior of £(u> A ) is similar to that represented in Figure 1 for the atom inside CN. Factor £(u> A ) is seen to rapidly 
decrease with raising the distance as it should be in view of the evident fact that photonic states coupled with CN 
electronic excitations (those increasing effective density p e ^{uj)) are spatially localized on the CN surface and their 
coupling strength with the excited atom decreases with raising the distance of the atom from CN. 

Figure || shows the ratio r r /r calculated according to Eq.(|^) for the atom in the centre of different CNs. (Note that 
r r /r = W s (u! A )/hu! A with W s (u! A ) being the total power of the atomic spontaneous radiation far from CN.) The ratio 
is very small indicating that the nonradiative decay dominates. However, the radiative decay is seen to essentially 
contribute in the vicinity of the interband transition frequencies so that the frequency dependence of W s {oj a ) (which, 
in principle, can be measured experimentally) reproduces CN electronic structure peculiarities. The main conclusion 
one can draw from Figure 3 is the Purcell effect in CNs, along with the increase of the atomic spontaneous decay rate, 
manifests itself by decreasing the power of the spontaneous radiation. 
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Our model of the atomic spontaneous decay in the presence of CN includes, as a limiting case, the model of the 
ideally conducting cylinder considered in Ref. 0. In particular, our Eq.(^J) reduces for the outward atomic position 
to Eqs.(15),(18) of Ref. [|| as a zz — * oo. The inset in Figure || shows factor £(uja) at ^ a = 3jo/h (kAR cn — 0.01) 
as a function of ta/Rcu for this case. The dependence is similar to that reported in Ref. |6j for z-oriented dipole at 
kARcn — 1- For the atom inside CN, Eq.(||) yields £(u>a) —> as a zz — ► oo — the result is natural since in this case 
only one eigen electromagnetic mode can propagate in CN; this mode is essentially transverse and, consequently, is 
not coupled with the atomic dipole moment oriented longitudinally. However, the actual £(wa) behavior discussed 
above is quite different from that predicted by the ideally conducting cylinder model since the latter one does not 
account for CN electronic quasiparticle excitations responsible for the nonradiative atomic decay dominating the total 
spontaneous decay process. 

Our theory may be generalized to cover the transverse atomic electric dipole orientation, electric quadrupole 
and magnetic dipole atomic transitions, properties of organic molecules inside/outside CNs |Q and of fullerene 
peapods p5| . The mechanism revealed of the photon vacuum renormalization is likely to manifest itself in other 
phenomena in CNs such as, for example, Casimir forces, electromagnetic fluctuations, etc. Jl6[ . 

The results of the present work may be tested by methods of atomic fluorescent spectroscopy and may possess 
various physical consequences. In particular, the effect of the drastic increase of the atomic spontaneous decay rate 
may turn out to be of practical importance in problems of the laser control of atomic motion |17| , yielding the increase 
of the ponderomotive force acting on the atom moving in the vicinity of CN in a laser field. One might expect the 
Purccll effect peculiarities that we predict for CNs to manifest in macroscopic anisotropically conducting waveguides 
with strong wave deceleration (for example, in microwave spiral or collar waveguides with highly excited Rydberg 
atoms inside). 

The work was financially supported by the Foundation for Basic Research of the National Academy of Sciences of 
the Republic of Belarus (Grants No F01-176 and No F02R-047). 
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Figure captions: 

Figure [jj: Factor £(u>a) calculated from Eq.(||) for the atom in the center of different "zigzag" CNs. Surface axial 
conductivities a zz appearing in Eq.(|^) were calculated in the r-approximation with t = 3 x 10~ 12 s _1 . 

Figure [2]: Factor £(w A ) for the atom at different distances outside "zigzag" (9,0) CN. Inset: £,{uja) at cl> j 4 = 37o/?i as 
a function of ta/Rch for the atom near (9,0) CN in the model of the ideally conducting cylinder. 

Figure S Ratio r r /T calculated from Eq.(@) for the atom in the center of different "zigzag" CNs. 
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